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PUNCHING TAPES: 



Analysis of Variance by Method of Fitting of Constants 
Complete program (SADOI or DOl) 

Approximate time in milliseconds 

Read-in: lOS^N 2 + N) 

Compute: 2(n - K)^ 

Output: 136 (N 2 + 3N + K 2 - K) 

where S = Number of samples 

N = Total number of variables 
K =' Number of dependent variables 
The program is read into the machine in the usual way. 
After, read-in a sum check is performed and if the tape 
has been read in correctly the machine will stop on the 
order 2k(0k2) l6 . Next the parameter tape is read in 
and the machine stops on 2k(oQk) l6 from 04k. The data 
tape is now read in and the machine prints out the 
required information for the analysis of variance; 
gross sums of squares and cross products, constants, 
inverse matrix and accounted for sums of squares and 
cross products . 

The total number of variables must satisfy the inequalities: 
N - K 

— — [3N - K + 3] + N < 792 and N < 30 

If K = 1, then N = 22. There is no limit on S unless the 
accumulated sums of squares exceeds 
Parameters needed are 

1. s = sample size. Punch sS on the tape. 

2. n = total number of variables. Punch nN on the tape. 

3. k = number of dependent variables. Punch kK on the tape. 
Immediately following the K there must be punched two 
carriage returns and then any identification desired, 
either letters preceded by a letter shift or numbers 
preceded by a number shift or any combination. At the 

end of the identification punch a 5th hole delay. 
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METHOD USED: 



Each observation is represented by a row of coefficients 
indicating the relationship of the observation to the 
constants to be fitted. These coefficients are each 
preceded by a + or - sign as appropriate, and the 
dependent variable s is the last item punched. These 
coefficients describe the complete classification of 
each observation; use a zero if the observation did 
not occur in a class and a 1 if it did. Constants for 
single degrees of freedom can be calculated by use of 
appropriate coefficients. Follow this by the terminating 
symbol N. If weights other than 1 are desired for any 
row of this observation matrix, precede the row by the 
numerical weight J punched as +JJ. The machine in effect 
treats these rows as if they occurred j times. 
If a terminating N is changed to F then the machine will 
stop so that another piece of the data tape can be 
inserted in the reader. 

The method is essentially that described in Kempthorne, 
"Design and Analysis of Experiments", page 79. The 
machine accumulates a sum of squares and cross products 
matrix from the observation matrix being read in. This 
matrix is inverted and solved for the various constants. 
The constants are printed out along with the inverse 
matrix for calculation of the standard errors of the 
constants . 

The standard errors of the constants can be calculated as 
follows. Subtract "Accounted for Sum of Squares" from 
"Raw Sum of Squares". This Residual Sum of Squares is 
* divided by S - (N - K) to provide Residual Variance. 
The diagonal elements of the inverse matrix multiplied 
by the residual variance give the Variances for each 
constant. The square roots of these Variances of Constants 
are the standard errors. 
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If k > 1, then the covariance model is in effect and the 
constants printed out are for the k variables in inverse 
order, i.e., the last variable first, next to the last 
2nd, etc. 

1. If the drum fails, the machine will stop on an FF 010 

from location 03K. 

ulated b ed 2^ th chin Hill 

stop on an FF(20l) 1 ^ from location 099- 
3. If any row of the observation matrix fails to have the 

correct number of variables the machine will stop on 

FFfeOO)^ from location 08F. 
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LOCATIOI I ORDER 



3* 
55 
36 
37 
38 

39 
40 

kz 

*5 

k6 
kj 
48 
*9 
30 



>*6 k6h 
k6 65L 
65L 
L3 (t)F 
L6 (x)F 
32 39L 
!*7 

50 7F 
L5 (x)P 
66 (t)F 
26 43L 
50 (x)F 
S3 F 
52 *2L 
50 110L 
75 (t)F 
66 (x)F 
47 46l 
4l 5F 
31 U5L 

ko zr 

L5 (x)F 

Ho r 

L5 (t)F 
40 F 
50 2F 
TJ IF 
L4 F 
40 (x)F 
L3 (x)F 
L6 5F 
56 51L 
L7 (x)F 
40 5F 



■OXEg 



PAG15 23 



if 'is" 



2^ f 
31 f 



27x55 
33S37 
12S53 
3^, I *2 T 



33,5* r 
28,5*' 



32* ,55 



24' Test size of loading element* 

31 

lo row interchange 
, Approximately zero in 

50» 
25 



29 1 



1-2' 



■39 



Bov interchange 
Address a parameter 



Store absolute 
value of largest 
element of row 



Linearly combine 
successive rows so 
as to get zero® 



UOCkttOS 


dUXB 


WOTB BU9B j 


51 


10 IF 








©almr the dlsgoml 




40 (t)F 


26,52 r 






of A 


58 


F5 51L 












^2 5H. 








■ ■ 


53 


42 45L 












L5 46% 








5* 


L4 109L 












40 48L 










55 


46 50L 












42 44l 










56 


L0 113L 












52 44L 










57 


L3 5F 




If s«ro doa f t ret c&li* | 




52 69L 








56 


LT. 5y 




> 1/2? 






32 61L 








59 


L5 66L 










46 65L 








6o 


26 64l 




J I* »(5) > 1/2 Multiply row 1/2 




F5 65L 




"I (scale dpva) 


61 


42 65L 










1.5 5F 








62 


00 IT 




*~ Determine if possible to ecajLe 19 




40 5F 




< 1/2 1(5) < X/k 


63 


LL 5F 










32 60L 








64 


50 7? 










L5 (x)F . 


66,28 f 






65 


10 (1)F 


3^,59' 








00 (1)F 


35,61 






66 


50 2F 




Waste 






40 (x)F 


67' -PQ 




Reacala rows 


67 


F5 66L 










42 66L 








68 


42 64L 










LO 114L 









LUUAXlUH 


UtU/CiK 




BOTES EMH&5 ! 


69 


56 6kL 




J ~ 1 




6f 








70 


6f 










L5 7? 






Determine If row 1 must haye further 


71 


L0 6F 






eliminations (dome 1 tlaes) 




36 2kL 






■ 


72 












ko jr 






I 


73 


LO 5P 






- Count number of rows I 




32 7^L 






■ 1 


7* 


22 8L 






Repeat for next row j 




kl 5F 






Set counter j 


75 


L5 1UJ. 






1 
1 




*0 (r)7 


6 




Set scaling factor i 


76 


L3 (y)F 






Terminate calcnlatiea if scaling < 2 ' j 




36 U1L 






p 

f 


77 


iH 6F 










U 7* 






Clear counters | 


78 


L5 *5L 






* E 




k2 88L 








79 


LO *F 






j 




Lfc 5F 








60 


1*2 83L 






| 




L5 88L 






- Set addresses I 


81 


LO kf 










k2 88L 


i 






82 


L5 83L 








- 


46 88L 




— 




83 


50 (y)F 










71 F 


8o,86 






8it 


ko F 










L5 85L 








89 


FO 7F 










LO kr 




















WW- 



LOGATIOH 


ORDER 


103 


4o (y-i)T 




F5 IF 


10k 


kO JF 




LO 3F 


105 


36 1019L 




22 80L 


106 


5F 




40 5? 


107 


LO k? 




32 108L 


108 


26 75L 




26 (iH2)F 


109 


00 IF 




00 IF 


110 


7L 4095? 




LL 4095F 


111 


kO F 




00 F 


112 


S6 21L 




41 (t)F 


113 


10 (t)F 




L3 F 


11* 


JO 2F 




40 (t)F 



90TSS 



«5 



95 



13 
14 



13' 



- Court a rovB 

Exit to Auxiliary II 
Repeat 



- Count m coltxH&s 

Repast 
Link 



•39 



1-2 



1/10 



■Sad constants 



LOCATION 



3 
k 

5 
6 

7 
8 

9 

10 
11 
12 

13 
14 

15 
16 

IT 



QKDH* 

06 1PQK 

K5 r F 
k2 QJL 
kl kF 
L5 10S5 
kO JF 
92 131F 
JO 37S9 
50 3L 
26 S9 
L5 0S5 
L4 10S5 
k2 12L 
42 72L 
kl 8F 
F5 10S5 
00 20F 
k6 11L 
46 53L 
50 S3 
50 9L 
26 S8 
00 SK 
00 F 
4l 3F 
92 131F 
L5 F 
40 OF ♦ 
F5 3F 
40 3F 
L5 OF 
00 IF 
56 13L 
L3 8F 
36 19L 
L5 3F 
L4 11S5 



BCTH3 



BASS 25 1 1 1$ 



Inverse - 
(H - K) + 1 



(N - K) + S 3 




To print out 




BAGS 29I 1: 1? 



TffiST&T 



0BS3SSR 



36 
37 

38 

39 
to 

*3 

*5 
to 

*7 
M 

49 
50 
51 
52 

53 



to TP 
10 *F 

LD 2S5 
32 9L 
4l IF 

F5 2S5 
to OF 
50 OF 
75 1065 
S5 F 
k2 IF 
L5 IF 

to 48L 
JO V3S9 
50 43L 
26 S9 
L5 1335 
to ^9L 
F5 to5 
to 5F 
kl 6F 
50 258J 
50A7L 
26 S8 

00 S7 
00 F 
15 89L 
to 52L 
41 3F 
50 S3 
50 51L 
26 S8 
00 SK 
00 F 
L5 88L 



(5 - K) + 1 



3k 
55 
56 
57 

58 

59 

60 

6l 

62 

63 
6k 

65 
66 

67 
68 

69 
70 
71 



571* 
kl k¥ 

kl IF 

50 IF 

00 IF 

S5 S3 

50 F 

7^ F 

IA IF 

1*0 IF 

LL IF 

36 68L 

U 9S5 

HO OF 

L5 56L 

1*2 62L 

1*2 63L 

L5 F 

10 UP 

1*0 F 

F5 63L 

1*2 63L 
1*2 62L 
F5 OF 
1*0 OF 
ID 10S5 
32 62L 
22 53L 
L5 51^ 
IA 5S5 
1*0 57L 
F5 kF 
kO kF 
LO 10S5 
32 56L 
35 OF 



LOCATIOK 


ORDER 


H0TES PACK 32 j 


72 


to OP 








L5 F 




JEL/ + | 


73 


to 2F 








F5 3F 




I 


7* 


to 3F 








L5 2F 




1 


75 


00 IF 




j 

1 




36 73L 






76 


92 131F 








92 515F 




j 


77 


L5 IF 








50 77L 




1 


78 


26 s6 




1 




F5 52L 






79 


F^ 10S5 

• 




i 




ik to5 




1 


80 


to 52L 








F5 6F 






81 


to 6f 








LO 5F 






82 


32 50L 








lA k&L 






83 


¥k 2S5 






to 48L 




84 


F5 5F 








to 5F 




1 


85 


92 131F 








92 515F 




i 


86 


L5 6F 




| 




LO 7S5 






87 


32 46L 








22 F 




Link | 


88 


30 S3 




i 




Jk 25S3 






89 


26 S8 








00 SK 1 







TBSXBM 


| ' QRB&R 


IOTE3 PASS 55 






00 F 
00 59F 

00 996K 




J »...-.* : , ^ 





L3 F 
3k 66F 




Sum check 


1 


FF F 
2k 66f 






2 


KB 2164F 
86 5859? 
26 L 
26 IN 








